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ABSTRACT 



Context. One of the main goals of the COROT mission is to get precise photometric observations of selected bright stars in 
order to allow the modelling of their interior through asteroseismology. However, in order to interpret the asteroseismological 
data, the effective temperature, surface gravity, and chemical composition of the stars must be known with sufficient accuracy. 
[ Aims. To carry out this task, we have developed a spectroscopic method called APASS (Atmospheric Parameters and 
■ Abundances from Synthetic Spectra) which allows precise analysis of stars with a moderate to high rotational velocity, which 
T-H is the case for most primary COROT targets. 

Methods. Our method is based on synthetic spectra in which individual lines are replaced by analysis units (isolated lines or 
line blends, depending on the crowding of the spectral region and on the rotational broadening). It works differentially with 
respect to the Sun and allows the atmospheric parameters and chemical abundances to be determined by considering analysis 
, units with different sensitivities to these various parameters. 
O I Results. Using high signai-to-noise spectra and the APASS method, we determined the atmospheric parameters and chemicai 
abundances of 13 primary COROT targets. Our resufts agree well with those obtained by Bruntt using his software VWA 
and with those obtained with the software TEMPLOGG. However, in both cases, our error bars are significantly smaller 
than those of other methods. Our effective temperatures are also in excellent agreement with those obtained with the IR 
, photometry method. For five stars with relatively low rotational velocity, we also performed an analysis with a classical 
equivalent -width method to test agreement with APASS results. We show that equivalent-width measurements by Gaussian 
or Voigt profile-fitting are sensitive to the rotational broadening, leading to systematic errors whenever the projected rotation 
velocity is non-negligible. The APASS method appears superior in all cases and should thus be preferred. 
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1. Introduction 

COROT (convection, Rotation and planetary Transits) 
will carry out several projects in asteroseismology 
(Baglin et al., 1998| ). It will also be the first satellite 
launched with the aim of detecting exoplanets by the tran- 
sit method ('Rou an et al., 2000| ). The instrument features 
a 27 cm telescope, and the asteroseismology program will 
use two 2048x2048 CCD cameras, leading to a field of 
view of 3.5 deg^. During the 2.5 years of the mission, 5 
fields of the Galactic plane will be observed continuously, 
each one during 150 days. Among these fields, a few bright 
stars have been selected as candidate primary targets of 
the mission, and getting a precise photometric monitoring 
of these stars to allow the modelling of their interior is the 
main goal of the asteroseismological part of COROT (see 
Table 1 for the preliminary list of these primary targets) . 

Send offprint requests to: M. Gillon 
email : gillon@astro.ulg.ac.be 



To allow thorough modelling of the stellar interior, as- 
teroseismological data are not enough, so additional con- 
straints must be provided by precisely determining the 
atmospheric parameters and abundances. To accomplish 
this task, we set up a spectroscopic analysis method, 
APASS (Atmospheric Parameters and Abundances from 
Synthetic Spectra), that aims at the highest possible pre- 
cision. With this method, we analysed the majority of the 
COROT main targets, including some stars later rejected 
from the preliminary main target list. 



In Sect. 2, we present the spectroscopic observations 
used in this work, and in Sect. 3, the method APASS is 
described. In Sect. 4, we present the results and compare 
them to previous studies. We give our conclusions in Sect. 
5. 
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Table 1. Preliminary list of COROT main targets from 
IBagiin et al., 2004| The V sini values were provided by C. 
Catala (private communication). HD 52265 has a giant 
planet detected by spectroscopy ( |Butler et al., 2000| ). 



Table 2. Observation log for the spectra of the stars we 
analysed. The signal-to-noise ratio (S/N) in the last col- 
umn is calculated around 6500 A. co-ad: co-addition of 
~ 100 spectra taken the same night. 



Star 


Spectral type 


V sini (km s ^) 




Star 


Instrument 


Date 


UT start 


S/N 


HD 52265 


GO HI-IV 


5.2 




HD 43587 


ELODIE 


14-Jan-98 


22:32 


250 


HD 49933 


F2 V 


10.9 




HD 57006 


ELODIE 


10-Dec-OO 


00:41 


250 


HD 49434 


Fl V 


89.6 




HD 45067 


ELODIE 


15-Jan-98 


22:20 


260 


HD 43587 


F9 V 


5.9 




HD 43318 


ELODIE 


19-Jan-98 


22:48 


120 


HD 170580 


B2 V 


16.0 


/3 Cep 


HD 49933 


ELODIE 


21-Jan-98 


22:17 


210 


HD 171834 


F3 V 


72.1 




HD 49933 


HARPS 


19-Jan-04 


co-ad 


750 


HD 177552 


Fl V 


41.2 




HD 55057 


ELODIE 


15-Jan-98 


23:10 


270 


HD 181555 


A5 V 


> 200 


S Scuti 


HD 184663 


ELODIE 


18-Jun-OO 


01:30 


170 


HD 180642 


B1.5 H-HI 




l3 Cep 


HD 49434 


ELODIE 


17-Jan-98 


22:33 


160 


HD 45067 


F8 V 


8.0 




HD 171834 


ELODIE 


05-Sep-98 


19:00 


300 


HD 46558 


FO V 


63.1 




HD 177552 


ELODIE 


14-Jul-98 


22:46 


260 


HD 43318 


F6 V 


7.5 




HD 175726 


ELODIE 


lO-Jul-01 


21:23 


140 


HD 175726 


G5 V 


13.5 




HD 46558 


ELODIE 


21-Dec-99 


01:08 


140 


HD 171234 


A3 V 


162.8 


5 Scuti 


HD 52265 


CORALIE 


24-Nov-Ol 


08:03 


140 



2. Spectroscopic observations 

We extracted a spectrum of each of the COROT main 
targets analysed here from the COROT ground-based 
asteroseismology database, Gaudi^ ( |Solano et al., 200'5| ). 
All but one were obtained with the ELODIE spectro- 
graph attached to the 1.93 m telescope at Observatoire 
de Haute-Provence (OHP). ELODIE is a fiber-fed cross- 
dispersed Echelle spectrograph that provides a complete 
spectral coverage of the 3800-6800 A region, with a resolv- 
ing power R = 45000 ( |Baranne et al., 1996| ). The spec- 
trum of HD 52265 was obtained with the CORALIE 
spectrograph attached to the 1.20 m Swiss telescope at 
La Silla (Chile), providing a resolving power R = 50000 
HQueloz et al., 2000| ). For one star, HD 49933, we also 
made an analysis with a co-addition of ^ 100 spectra 
obtained using the ultra-high resolution (R = 115000) 
HARPS spectrograph attached to the 3.60 m ESO tele- 
scope at La Silla ( [Mayor et al., 2003| ). Table 2 gives the 
observation log for the spectra used in this analysis. 



3. Description of the method 

Using the assumptions of Local Thermodynamic 
Equilibrium (LTE) and plane-parallel, homogeneous 
atmosphere, APASS is a differential method based on the 
synthesis of selected spectral regions with Sun-calibrated 
line data and comparison with a stellar spectrum in 
order to determine the projected rotational velocity (v 
sini), the effective temperature Teg, the surface gravity 
log g, the microturbulence velocity eturb, as well as the 
atmospheric abundances. 



^ http://sdc.laeff.esa.es/gaudi 



3.1. Analysis Units (Anil) 

We developed APASS in order to be able to precisely 
analyse stellar spectra with a moderate to high projected 
rotational velocity (see Table ^ , which cannot be anal- 
ysed with a classical equivalent -width (EW) method. Such 
spectra are convolved with the rotation profile, and most 
lines are blended with other ones, forming "packs" of lines. 
One of the main characteristic of APASS is that it con- 
siders either individual lines or "packs" of lines, treating 
them as analysis units (AnU). Those AnU are selected in 
the first step of the analysis, by visual inspection of the 
stellar spectrum. All the AnU that are dominated by iron 
(iron-dominated AnU : Fe-AnU) are used in determining 
the atmospheric parameters, while the other AnU are used 
to determine the atmospheric abundances. 

3.2. Atmospheric models 

We used models interpolated in the Vienna model 
atmosphere grid ( [Nendwich et al., 20031 120041 
[Heiter et al., 2002 ) . They were computed with a modified 
version of the ATLAS9 code ( |Kurucz, 1993| ) in which 
turbulent convection theory in the CGM formulation 
( Canuto et al., 1996| ) is implemented. 

3.3. Atomic data 

APASS uses limited spectral regions selected on the basis 
of the following criteria : 

— scarcity of telluric lines, 

— absence or scarcity of molecular lines, 

— high diversity of atomic lines: different elements, vari- 
ous sensitivities to the atmospheric parameters. 
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— presence of numerous pseudo-continuum windows to 
accurately normalize the continuum in order to gain 
precision. 

The atomic lines data were extracted from the VALD 
database (|Kupka et al., 1999| [Piskunov et al., 1995| 
IRyabchikova et al., 1999|) and the Kurucz atomic lines 
database ^Kuruc z*, 1998). For all the lines with a 
significant equivalent width, we made a new deter- 
mination of the log gf from an inverted analysis of 
the Jungfraujoch Solar Atlas | |Delbouille et al., 1973| ), 
using a CGM solar model having {T^s, log g, eturb) = 
(5778 K, 4.44 dex, 1.00 km s~ ^) and solar abundances 
from Grevesse & Sauval (1998) Whenever possible, we 



used damping constants calculated following Anstee & 
1199511 ■ which have proved to be quite accurate 



O'Mara 

HAnstee et al., 1997| [Barklem fc O'Mara, 2000| ). 

3.4. Parameter determination 

APASS is an iterative method. Starting from an initial 
model, the method consists of cycles of determination of 
the parameters, repeated until convergence for all param- 
eters. Each cycle includes adjusting the continuum, de- 
termining the level of convolution of the spectrum, and 
determining the atmospheric parameters Teffj log g, eturb, 
and [Fe/H]. The projected rotational velocity v smi is de- 
termined afterwards. 

3.4.1. Initial model 

For each star analysed in this work, we begin with an ini- 
tial atmospheric model computed with parameters T^g, 
log g, eturb and a metallicity extracted from the litera- 
ture or obtained using Stromgren indices from the Gaudi 
database and the sofware TEMPLOGG ( [Rogers, 1995| ). 
To synthesize the spectrum of the selected regions, we also 
need to set initial values for two parameters representing 
the convolution of the spectrum: 

— the projected rotational velocity v sim. As initial value, 
we use the v sini determined by C. Catala (private 
communication) , 

~ the parameter representing the combined effect of 
the instrumental profile and the stellar macroturbu- 
lence. We use an initial value of 4 km s~^ for fl and 
assume a Gaussian shape for that convolution profile. 



3.4.2. Correction of the Doppler shift 

Before starting the first cycle of parameter determination, 
we use the synthetic spectrum computed with the initial 
model to determine the Doppler shift of the stellar spec- 
trum. For this purpose, we select a large number of AnU 
and determine the ratio of their position in the stellar 
spectrum and in the synthetic spectrum. We use the av- 
erage ratio for correcting the spectra that are synthesized 
afterwards. 




0. 2 - 

5321 



5322 
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Fig. 1. Portion of the stellar spectrum of HD 52265 
(empty circles) showing two Fe-AnU. Overplotted: syn- 
thetic spectrum computed with ~ 0.75 and 6.75 km 
s~^ (dotted lines) and with = 3.75 km s^^, the best-fit 
value (continuous line), lambda = wavelength in A, flux 
— normalized stellar flux. 



3.4.3. Adjustment of the continuum 

We begin each cycle with a new determination of the con- 
tinuum of the observed spectrum on the basis of a number 
of pseudo-continuum windows selected from inspection of 
both the observed and synthetic spectra. The mean flux 
in these windows is measured in both the theoretical and 
the observed spectra, and a table is constructed that con- 
tains the ratios of the mean fluxes. A spline curve is then 
fitted through these points, and the stellar spectrum is 
divided by that curve, providing our re-normalized spec- 
trum. This method allows us to take into account the weak 
lines present in the pseudo-continuum windows as accu- 
rately as possible. 

3.4.4. Determination of the convolution parameter 

The convolution parameter Q (macroturbulence + instru- 
mental profile) is determined for every Fe-AnU by ad- 
justing the synthetic spectrum with an minimization 
method, and the average overall Fe-AnU is the adopted 
value for fl (see Figure^. During all iteration cycles, we 
keep V sini fixed at its initial value. Indeed, as these two 
broadening parameters appear to be strongly anticorre- 
lated in an minimization, allowing both of them to 
vary produces instabilities in the algorithm. 

3.4.5. Determination of the Fe abundance 

A value of the Fe abundance is deduced from each Fe-AnU 
in the following way. First, the user defines the limits in 
wavelength of the Fe-AnU considered. Then, the surface 
under the Fe-AnU is computed in between the two lim- 
its of the observed spectrum. This is a kind of "observed 
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EW" . The same is done for the synthetic spectrum. If the 
two surfaces are not equal, the Fe abundance is modified 
and a new synthetic spectrum is computed, and so on un- 
til the computed surface matches the observed one. The 
adopted Fe abundance is the average of the values found 
from all the Fe-AnU considered. 

3.4.6. Determination of the atmospheric parameters 

eturb, Tcff, and log g 

In a first step, the sensitivity of each Fe-AnU to the atmo- 
spheric parameter considered is determined. First, we deal 
with the microturbulence velocity eturb- We adopt a value 
of eturb a-iid compute the Fe abundance from each Fe-AnU. 
Then, we change the value of eturb by a fixed amount and 
recompute the Fe abundances using this new model. The 
difference in Fe abundance obtained when changing eturb 
gives the sensitivity of each Fe-AnU to this parameter. 
Once the sensitivity of each Fe-AnU to eturb is determined, 
a diagram of abundance versus sensitivity is constructed 
and eturb is determined in order to remove any correlation: 
the deduced abundance should not depend on the sensi- 
tivity to eturb- This procedure is similar to the classical 
method where eturb is adjusted so that the abundance does 
not correlate with the line EW (the latter roughly giving 
the sensitivity to eturb for weak to medium strength lines). 
The microturbulence is determined in order to remove any 
correlation between the abundances and the sensitivities 
(see Figure El . In simple words, we want all the Fe-AnU 
to give us the same global abundance, whatever their sen- 
sitivity to eturb - 

Following the same method, the correction to the 
abundances of our Fe-AnU is computed for two models 
differing by their value of Toff. Once the sensitivity of our 
Fe-AnU to Tgg is computed, T^g is adjusted to remove 
any correlation between the abundances and the sensitiv- 
ities to Toff, which is similar to removing the trend in the 
abundance versus line excitation potential plot. 

The same method is applied again to determine log g; 
i.e. the sensitivity to log g of the Fe-AnU is computed, and 
the surface gravity is adjusted to remove any correlation 
between the abundances and the sensitivities to log g. 

A complete cycle includes determinations of fl, [Fe/H], 
eturb, Toff, and log g. The iteration ends after convergence, 
i.e. when the corrections on all parameters become much 
smaller than the error bars. 

3.4.7. Determination of v 5\ni 

Once convergence is reached, v sini is determined by ad- 
justing the synthetic spectrum by an minimization for 
every Fe-AnU, and the average value is taken as the fi- 
nal V sini. If this final value significantly differs from the 
initial value, the whole iteration process is repeated. As 
mentioned above, v sini is not adjusted during the iter- 
ation cycles because it is strongly anticorrelated with fl. 
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Fig. 2. Plot of the abundance deduced from each AnU as 
a function of the sensitivity of the AnU to eturb , for eturb 
— 1.20 km s^^ (top) and for eturb = 1-30 km s^^ {bottom). 
The slopes obtained by linear regression are -0.96 ± 0.24 
for eturb = 1-20 km s~^ (initial guess) and -0.04 ± 0.24 for 
eturb = 1-30 km s~^ (adopted value). 

Adjusting both parameters together thus produces insta- 
bilities. 

3.5. Computation of the error bars on the parameters 

For each parameter, we must distinguish the intrinsic un- 
certainty and the uncertainties coming from errors on the 
other parameters (external uncertainties). 

3.5.1. Intrinsic uncertainties 

For Toff, log g, and eturb, the standard deviation on the 
slope in the sensitivity versus abundance of Fe-AnU di- 
agram allows us to compute the intrinsic uncertainty on 
these parameters, i.e. the statistical deviation. As an ex- 
ample, the uncertainty on the slopes in Figure El is 0.24, 
and the difference between the two slopes is 1.00, for a dif- 
ference in eturb = 0.10 km s~^ between the models. The 
intrinsic uncertainty on eturb will be in this case: 

CT = 0.10/1.00 X 0.24 = 0.02 km s-i. (1) 
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The convolution parameters v sini and 17, as well as 
the iron abundance [Fe/H], are determined independently 
for each Fe-AnlJ. The adopted value is the average value 
for all the Fe-AnU used, and the intrinsic uncertainty is 
the standard deviation of the mean. 

3.5.2. Uncertainty on one parameter due to the 
intrinsic uncertainty on another parameter 

To compute the uncertainties on, e.g., log g coming from 
another parameter, e.g.. Toff, we make a new determina- 
tion of log g using a model with Tcff differing from the 
adopted value by its intrinsic uncertainty. As an exam- 
ple, if we have a final value for log g — 4.34 dex and a 
final value for Tcs = 6400 K with a standard deviation — 
50 K, we set Toff to 6450 K (or 6350 K) and determine 
log g. The difference between the obtained value for log g 
and 4.34 dex is the uncertainty on log g coming from the 
uncertainty on Teg. 

For each atmospheric parameter, we compute the un- 
certainties coming from all the other parameters. The to- 
tal error bar is then computed as a quadratic sum of these 
uncertainties. The intrinsic uncertainty is added quadrati- 
cally to the external uncertainties to obtain the final error 
bar. An exception exists for v sini, for which we consider 
the intrinsic uncertainty and the uncertainty coming from 
17 as fully correlated. In this case, the intrinsic uncertainty 
is added linearly to the uncertainty coming from O, and 
the result is added quadratically to the other uncertainties 
to obtain the final error bar on v sim. 

3.6. Test of the APASS method on an ELODIE 
spectrum of the Moon 

As a test, we determined the parameters of the Sun us- 
ing an ELODIE spectrum of the sunlight reflected by 
the Moon. Starting with T^s = 6000 K, log g = 4.00 , 
eturb = 2.00, [Fe/H] = 1.00, and v sim = 5 km s~^ as 
initial guess parameters, the method converged to the fol- 
lowing parameters: Tcff = 5786 ±25 K, log g = 4.45 ±0.03, 
et^,,b = 1.01 ± 0.06 km s^^ , [Fe/H] = +0.01 ± 0.02, and v 
sini — 1.94 ± 0.72 km s^^. These are very close (and well 
within error bars) to the "expected" solution, i.e. the set of 
values used in our solar model: Toff = 5778 K, log g = 4.44, 
et^^.b = 1-00 km s-\ [Fe/H] = 0.00 , v sim = 2.00 km a'^. 
This test demonstrates the correct convergence of APASS 
and the consistency of the ELODIE spectrum with the 
Jungraujoch Solar Atlas. 

3.7. Abundance determination 

Once we have obtained the atmospheric parameters, we 
determine the atmospheric abundances for every element 
for which we have AnU containing significant lines of this 
element. We proceed by steps, in order to avoid accumu- 
lating the uncertainties. We first adjust the abundance of 
elements within AnU containing lines of these elements 



alone or blended only with iron lines. For the sake of clar- 
ity, let us call A such an element. Once we have determined 
the Fe abundance and the abundance of A, we can pro- 
ceed determining the abundance of an element X whose 
lines are blended with lines of Fe and A. And so on, until 
all elements with significant lines are included. 

3.8. Computation of the error bars on the abundances 

The uncertainty is computed taking into account the part 
coming from the uncertainty on the abundances of the 
other elements included in the AnU that we used, along 
with the uncertainty on the parameters. Let us consider an 
example. We determine [Si/H] using ten lines, contained 
in 8 different AnU. The standard deviation of the mean 
is 0.02 dex. Eight of these lines are blended in their AnU 
with Fe and Cr lines. We quadratically sum 0.02 dex to 
the intrinsic dispersion on [Cr/H] and [Fe/H], obtaining 
e.g. 0.04 dex. Then with different atmospheric models we 
determine the uncertainties on [Si/H] coming from the in- 
trinsic uncertainties on the atmospheric parameters; i.e. 
we redetermine [Si/H] with a model whose T(.g equal the 
adopted value plus or minus the intrinsic uncertainty on 
Toff, and we do the same for log g, eturb, 17, and v sini. 
We then quadratically sum all the uncertainties on [Si/H], 
obtaining the final error bar. For elements for which we 
do not have enough AnU to get a significant standard de- 
viation of the mean, we take the average of standard devi- 
ations for the other elements except iron. We set to 4 the 
minimum number of AnU needed to judge the standard 
deviation significant. 

4. Results - comparison 

4.1. Stellar parameters 

The final stellar parameters are presented in Table 4. The 
results are compared to other determinations: 

- We used the software TEMPLLOG ( [Rogers, 19951 ) 
and the Stromgren indices extracted from the Gaudi 
database (Table 3) to determine the atmospheric pa- 
rameters Tcff, log g, and [Fc/H], with a precision 
around 200-250 K, 0.3 dex, and 0.2 dex, respectively, 
according to Rogers. For HD 52265, the H^ index was 
not available on Gaudi, so we took it from Hauck & 
Mermilliod (1998). 

— For nine of the stars, Tcs was determined with the IR 
photometry method developed by Ribas et al. H2003|l . 
considered as more reliable than TEMPLOGG by the 
COROT team (C. Catala, private communication). 

- For the stars HD 43587, HD 43318, HD 45067, HD 
57006, and HD 49933, which have a low enough ro- 
tation velocity, we also used our equivalent-width 
method to analyse them. This method is described in 
Bruntt et al. (2004). 

— The star HD 52265 has a massive planet orbiting 
it that was detected by the radial velocity method 
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Table 3. Stromgren photometric indices of the COROT 
main targets analysed in this work taken from the GAUDI 
database. The index for HD 52265 is taken from Hauck 
& MermiUiod (1998) 



Star 


V 


H/3 


(b-y) 


ml 


cl 


HD 52265 


6.294 


2.619 


0.360 


0.186 


0.423 


HD 43587 


5.720 


2.633 


0.382 


0.188 


0.366 


HD 43318 


5.629 


2.656 


0.320 


0.154 


0.462 


HD 45067 


5.880 


2.638 


0.364 


0.165 


0.405 


HD 57006 


5.916 


2.634 


0.336 


0.181 


0.494 


HD 49933 


5.782 


2.659 


0.268 


0.133 


0.466 


HD 175726 


6.731 


2.621 


0.366 


0.176 


0.312 


HD 177552 


6.510 


2.675 


0.243 


0.132 


0.545 


HD 184663 


6.430 


2.675 


0.278 


0.149 


0.466 


HD 46558 


6.887 


2.736 


0.244 


0.143 


0.542 


HD 171834 


5.416 


2.675 


0.250 


0.148 


0.543 


HD 49434 


5.753 


2.759 


0.170 


0.181 


0.705 


HD 55057 


5.463 


2.742 


0.184 


0.182 


0.877 



jButler et al., 2000| ). Santos et al. derived stellar pa- 
rameters for planet-host stars using an equivalent- 
width method l |Santos et al., 2004| ). For HD 52265, 
they made two analyses, one based on a CORALIE 
spectrum and the other on a FEROS spectrum. 
— Nine of the stars analysed in this work were 
also analysed by Bruntt using the VWA software 
HBruntt et al., 2002||Eruntt et al., 200^ . This method 
is based, as is APASS, on the comparison of the stellar 
spectrum with synthetic spectra. 

Figure 3 shows the comparison of our values for Toff, log 
g, and [Fe/H] with the values obtained with TEMPLOGG 
and our equivalent-width method, and with the values 
obtained by Bruntt with VWA, while Figure 4 shows the 
comparison of our values for Toff with the values obtained 
with the IR photometry method. 

4.1.1. Comparison with TEMPLOGG values 

For Toff, we found an average difference < S > = -2 
K when comparing our values with those obtained with 
TEMPLOGG. Nevertheless, as can be seen in Figure 3, 
the dispersion is high {a = 253 K). Despite the small 
number of points, a tendency seems to appear in the plot. 
TEMPLOGG tends to give higher values than APASS for 
Toff ~ 6000 K, and lower values for higher Toff- For log g, 
<5>= -0.03 dex with a ^ 0.30 dex, and for [Fe/H] the 
< S >= +0.07 dex and a — 0.18 dex. We thus conclude 
that our results generally agree with the values obtained 
with TEMPLOGG. However, the scatter is rather large, 
and can generally be explained by the larger error bars in 
the photometric TEMPLOGG determination. 

4.1.2. Comparison with IR photometry values 

As can be seen from Figure 4, agreement between the Toff 
obtained with APASS and from IR photometry is very 
good, within the small error bars for most of the stars: 




6000 6500 7000 

Teff(K) APASS 



Fig. 4. Comparison of the Tcfr derived in this work with 
those obtained with the IR photometry method of Ribas 
et al. (2003). 

< 5 >= +40 K, with a = 110 K. Only the hottest star, 
HD 49434, shows significant disagreement. If we do not 
consider it, < 5 > is reduced to -1-9 K and a to 55 K. 

4.1.3. Comparison with the values obtained with our 
EW method 

We find < S >= +167 K with cr = 187 K for Teff, 
when comparing our results using APASS with those ob- 
tained using our EW method. The values obtained for 
[Fe/H] (< 5 >= -f0.09 dex, a = 0.11 dex) and for etmb 
(< 6 >= —0.75 km/s, a — 0.86 km/s) are not in very 
good agreement either. Only the values obtained for log g 
are consistent (< S >= 0.00 dex, a = 0.08 dex). 

Our EW method uses Gaussian and Voigt profile- 
fitting to measure the EW of single lines. All the stars 
analysed here have a v sini higher than the Sun, and we 
wanted to know if the measured EW are independent of 
the rotational broadening. We measured the EW of three 
synthetic Fel lines computed with different convolution 
parameters v sim, with solar values for eturb, T^g, log g, 
and [Fe/H] and with a fixed value of fl. We also computed 
the theorical equivalent-width of these lines. Our results 
in Figure 5 show that the measured EW is sensitive to 
the rotational broadening of the lines. The fit of Gaussian 
profiles tends to give EW values that are too low for low v 
sini, and the measured EW tends to increase with v sini. 
The fit of Voigt profiles gave good results for low values of 
V sini, but seemed to give absolutely incorrect values for 
values of v sini higher than 4 km/s. The same test with a 
fixed value of v sini and different values of Q showed also a 
clear dependence of the measured EW on the convolution 
magnitude of the spectrum. 



M. Gillon & P. Magain: Atmospheric parameters of COROT main targets 



7 




6000 6500 








-0.2 



— -0.4 





111111/ 




/ ++ - 






- / 4 

/] 1 


1 1 1 1 1 



-0.4 -0.2 



< 7500 
^ 7000 
H 6500 
^ 6000 




J I I 

6000 7000 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
3 3.5 4 4.5 5 



> 




o 8000 

o 

S 7500 



0) 



7000 
6500 
6000 




6000 8000 
Teff(K) APASS 




o o.b 



3.5 4 4.5 5 
log g APASS 



o 
a 
o 



0) 




-0.5 0.5 
[Fe/H] APASS 



Fig. 3. Comparison of the atmospheric parameters Toff, log g , and [Fe/H], derived in this work with those obtained 
with our EW method {top), by Bruntt using VWA {middle), and with TEMPLOGG using Stromgren indices from 
the Gaudi database {bottom). 



We plot in Figure 6 the observed differences between 
the two methods as a function of v sim. For T^g, a clear 
tendency appears. The higher the rotational broadening 
of the spectrum, the larger the discrepancy between the 2 
methods, with the EW method giving values lower than 
APASS. For the other parameters, we can only say that 
lower [Fe/H] and higher eturb are found with the EW 
method. All those results seem to indicate that an EW 
method using Gaussian and Voigt profiles fitting is not 



suitable for the analysis of stars with a projected rotation 
velocity that is significantly different from that of the Sun. 

4.1.4. Comparison with VWA values 

A small < ,5 > = + 29 K is found for T^g (o" = 140 K), 
and < S > = -0.02 dex for log g with a. a = 0.21 dex. 
For [Fc/H], there is marginal disagreement, with < S > = 
+0.10 dex and a = 0.13 dex. Despite the fact that these 
two methods use synthetic spectra to obtain atmospheric 
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Fig. 5. Measured equivalent -widths of 3 synthetic Fel hues by Gaussian (sohd hue) and Voigt (dashed hue) profile- 
fitting, for different values of v sini. is fixed to 2 km/s. The theorical value of the EW appears in the upper left 
corner of each plot. 



parameters, our error bars are significantly smaller than 
those obtained with VWA. We believe that this could be 
due to four causes : 

1. strict selection of the spectral regions used for the anal- 
ysis, 

2. precise determination of the continuum, 

3. precise determination of the convolution level of the 
spectrum, 

4. the use of packs of blended lines as analysis units, 
with all the log gf of the main lines of the packs pre- 
cisely determined from a high resolution solar atlas. 
This point is very important. If we would restrict our 
analysis to isolated lines, the number of independant 
data available for our analysis would decrease drasti- 
cally with the increase of the line broadening. Using 
blended lines allows us to extract a high quantity of 
information from spectra, even if they are moderately 
broadened. 



4.1.5. HD 52265: comparison with results obtained by 
Santos et al.(2004) 

Table 4 shows that our results for HD 52265 with APASS 
are in good agreement with those of Santos et al. from the 
analysis of their FEROS spectrum of this star, as can be 
seen from Table 4. As our analysis is based on a CORALIE 
spectrum, it is more meaningful to compare our results 
with the CORALIE results of Santos et al. We find ATes 
= 103 K. As can be seen from Figure 6, it agrees with the 
overall tendency noticed when comparing the results of 
APASS and of our EW method. A linear regression made 
on the values of v sini as a function of the difference in 
Teff gives an intercept = 1.54 km s^^, near the observed 
V Sim of the Sun, 1.9 ± 0.3 km s'^ ( jSoderblom, 1982| ). 
The correlation coefficient r is 0.97. Our values for eturb, 
log g, and [Fe/H] are within the error bars of Santos et 
al., but their values are not inside our own error bars. For 
these parameters, no information comes from Figure 6, 
except the fact that the EW method used by Santos et 
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Fig. 6. Differences between the values obtained with APASS and with the EW method as a function of v sini. Open 
symbol: results obtained by Santos et al. for HD 52265 with a CORALIE spectrum. 



al. docs not seem to give much higher eturb than APASS, 
unlike our own EW method. We notice that Santos et al. 
measured the equivalent-width by fitting a Gaussian pro- 
file, whatever the strength of the line. This method will 
tend to give too low an EW for strong lines, while Cturb 
is determined in the EW method by cancelling the cor- 
relation between the line abundances and their EW. An 
EW that is too low for stronger lines will produce an un- 
derestimation of eturb, so the fact that the eturb found by 
Santos et al. for HD 52265 is near our own value could be 
a coincidence that comes from the combined effects of {i) 
the systematic error due to the use of Gaussian profiles, 
whatever the strength of the line, and of (ii) the system- 
atic error due to the fact that the measurement of the EW 
by non-rotationally broadened profiles fitting is sensitive 
to the degree to which the spectrum is convolved. 



4.2. Abundances 

We present the abundances obtained in this work in Table 

7. These are abundances relative to the Sun. From the 
abundance patterns we find no evidence of chemically pe- 
culiar stars. 



5. Conclusions 

We have developed a method able to precisely determine 
the atmospheric parameters and abundances of stars with 
moderate to high projected rotational velocity. This dif- 
ferential method, APASS, is based on the plane-parallel 
homogeneous atmosphere and LTE hypothesis. 

With APASS we performed a parameter determina- 
tion and a detailed abundance analysis of thirteen poten- 
tial COROT main targets. The precision of our results is 
very high. Four stars in the preliminary list of COROT 
main targets could not be analysed with APASS. The S 
Scuti stars HD 181555 and HD 171234 have a projected 
rotational velocity that is too high for our spectroscopic 
method. For v sini larger than 150 km/s, it becomes very 
difficult to adjust the continuum correctly and to find a 
significant number of well-defined AnU, making APASS 
unsuitable for spectroscopic analysis. On the other hand, 
the (3 Cep stars HD 170580 and HD 180642 arc too differ- 
ent from the Sun for our differential method. Furthermore, 
LTE becomes a poor assumption for such hot stars. 

Prom the abundance pattern of the analysed stars, 

we found no evidence of chemical peculiarities. For most 
of the stars, our results for the fundamental parameters 
agree both with the estimates from TEMPLOGG us- 
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ing Stromgren photometry and with results obtained by 
Bruntt using his software VWA. The precision of our re- 
sults is high, although we take the influence of the uncer- 
tainty on the broadening parameters into account in our 
computation of the error bars. 

We also made an analysis of five COROT main targets, 
each with a moderate rotational velocity, using a classical 
equivalent-width method. The results show a clear dis- 
crepancy between the two methods. We than investigated 
this discrepancy and found that a classical equivalent- 
width method using Gaussian or Voigt profile-fitting is 
not adapted to the analysis of stars with a projected rota- 
tional velocity that is significantly different from the solar 
one. 

For the star HD 52265, which has a massive planet 
orbiting around it and which was analysed by Santos et 
al. (2004) using a classical equivalent-width method, our 
results agree well with the results obtained by these au- 
thors. Nevertheless, comparison of the results obtained af- 
ter analysing spectra coming from the same instrument 
(CORALIE) seems to confirm that an equivalent-width 
method using Gaussian or Voigt profiles tends to give sys- 
tematic errors on the parameters and abundances due to 
the sensitivity of the measured equivalent-width to the 
broadening of the spectrum. The use of rotationally broad- 
ened profile could lead to more reliable results, but only 
in the case of a rather slow rotation. For stars with a 
significantly higher projected rotational velocity than the 
Sun's, the concept of isolated lines is lost and only syn- 
thetic spectra methods such as APASS should be used 
to spectroscopically determine the atmosphere parameters 
and abundances. 
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Table 4. COROT targets and stellar parameters obtained with APASS. HD 49933 was analysed with an ELODIE 
spectrum (E) and a high S/N HARPS spectrum (H). See text for more details and Table 5 for references. 



Star 


V sini (km/s) 


eturb (km s"^) 


Teff (K) 


log g (dex) 


[Fe/H] (dex) 


Ref. 


HD 52265 


5.0 ±1.3 


1.30 ±0.04 


6179 ± 18 


4.36 ±0.03 


+0.24 ±0.02 


1 






1.33 ±0.08 


6131 ±47 


4.35 ±0.13 


+0.25 ±0.06 


3 






1.38 ±0.09 


6076 ± 57 


4.20 ±0.17 


+0.20 + 0.07 


4 








6135 ± 200 


4.09 ± 0.20 


+0.05 ± 0.20 


5 








6136 ± 51 






7 


HD 43587 


5.8 ±0.9 


1.18 ±0.04 


5978 ± 18 


4.29 ±0.02 


+0.02 + 0.02 


1 








5870 ± 60 


4.20 ±0.15 


-0.09 + 0.11 


6 






2.13 ±0.10 


5867 ± 20 


4.29 ± 0.06 


-0.13 ±0.02 


2a 








6327 ± 200 


4.54 ± 0.20 


+0.24 ± 0.20 


5 








5925 ± 65 






7 


HD 43318 


7.6 ± 2.2 


1.56 ± 0.08 


6196 ± 33 


3.68 ± 0.06 


-0.16 ±0.03 


1 








6190 ± 90 


3.70 ±0.15 


-0.21 ±0.11 


6 






2.13 ±0.31 


6045 ± 40 


3.63 ± 0.07 


-0.24 ± 0.03 


2a 








6517 ±200 


4.25 ± 0.20 


-0.06 ± 0.20 


5 








6220 ± 98 






7 


HD 45067 


7.9 ± 1.3 


1.31 ±0.05 


6070 ± 23 


3.99 ± 0.04 


-0.02 ±0.02 


1 








5970 ± 100 


3.80 ± 0.15 


-0.17 ±0.11 


6 






2.00 ±0.13 


5913 ± 30 


3.91 ± 0.08 


-0.15 ±0.02 


2a 








6358 ± 200 


4.37 ± 0.20 


+0.02 ± 0.20 


5 








6092 ± 59 






7 


HD 57006 


8.2 ± 1.5 


1.68 ± 0.06 


6177 ±26 


3.66 ± 0.05 


-0.02 ±0.02 


1 








6180 ± 70 


3.60 ±0.15 


-0.08 ±0.11 


6 






2.72 ± 0.09 


6026 ± 22 


3.61 ± 0.05 


-0.09 ± 0.02 


2a 








6250 ± 200 


3.74 ±0.20 


+0.06 + 0.20 


5 


HD 49933 


10.3 ± 1.8 


1.65 ±0.10 


6700 ± 65 


4.25 ±0.09 


-0.34 + 0.03 


1 (E) 




10.0 ± 1.6 


1.84 ±0.08 


6735 ± 53 


4.26 ±0.08 


-0.37 + 0.03 


1 (H) 








6780 ± 70 


4.30 ± 0.20 


-0.30 + 0.11 


6 






2.53 ±0.31 


6467 ± 56 


4.27 ± 0.08 


-0.37 ± 0.04 


2a (E) 






2.12 ±0.24 


6539 ± 47 


4.41 ± 0.07 


-0.45 ±0.02 


2a (H) 






2.33 ± 0.24 


6538 ± 44 


4.38 ± 0.08 


-0.47 ± 0.02 


2b (H) 








6543 ± 200 


4.24 ± 0.20 


-0.38 ± 0.20 


5 








6747 ± 64 






7 


HD 175726 


12.7 ± 1.7 


1.32 ±0.08 


6217 ±32 


4.61 ± 0.04 


+0.03 ± 0.03 


1 








6227 ± 200 


4.73 ± 0.20 


-0.05 ±0.20 


5 


HD 177552 


38.9 ±2.5 


1.88 ±0.10 


7025 ± 55 


4.45 ± 0.08 


-0.03 ± 0.03 


1 








6657 ± 200 


4.03 ± 0.20 


-0.38 ± 0.20 


5 








7042 ± 69 






7 


HD 184663 


57.3 ±5.3 


2.01 ±0.26 


6511 ± 108 


4.28 ± 0.21 


-0.13 ±0.06 


1 








6600 ± 90 


4.50 ±0.50 


-0.19 ±0.21 


6 








6700 ± 200 


4.40 ± 0.20 


-0.13 ±0.20 


5 


HD 46558 


60.0 ± 5.9 


1.63 ±0.21 


7053 ± 92 


4.44 ±0.14 


+0.04 ± 0.07 


1 








7223 ± 200 


4.71 ± 0.20 


-0.24 ± 0.20 


5 








6922 ± 66 






7 


HD 171834 


70.0 ± 6.3 


1.64 ±0.21 


6833 ± 105 


4.14 ±0.15 


+0.01 + 0.07 


1 








6840 ± 200 


4.60 ± 0.50 


-0.25 + 0.21 


6 








6658 ± 200 


4.04 ± 0.20 


-0.17 ±0.20 


5 








6879 ± 67 






7 


HD 49434 


85.4 ± 6.6 


2.77 ±0.26 


7632 ± 126 


4.43 ± 0.20 


+0.09 + 0.07 


1 








7300 ± 200 


4.40 ± 0.45 


-0.04 + 0.21 


6 








7343 ± 200 


4.21 ±0.20 


+0.01+0.20 


5 








7346 ± 69 






7 


HD 55057 


130.1 ± 7.5 


3.20 ± 0.27 


7478 ± 105 


3.84 ±0.17 


+0.26 + 0.08 


1 








7580 ± 250 


3.60 ± 0.50 


+0.14 + 0.21 


6 








7141 ± 200 


3.47 ±0.20 


+0.08 ± 0.20 


5 
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Table 5. The references for Table 4 



Ref.l This study with APASS 

Ref. 2a EW method, Kurucz 1993 model 

Ref. 2b EW method, CGM model 

Ref. 3 Santos et al. (2004), EROS spectrum 

Ref. 4 Santos et al. (2004), CORALIE spectrum 

Ref. 5 TEMPLOGG using Stromgren indices from Gaudi database 

Ref. 6 Bruntt et al. (2004) using VWA 

Ref. 7 IR photometry (Ribas et al., 2003) 



Table 6. Abundances obtained with APASS for the COROT main targets. HD 49933 was analysed with an ELODIE 
spectrum (E) and a high S/N HARPS spectrum (H). See text for more details. 



HD 52265 HD 43587 HD 43318 HD 45067 HD 57006 HD 49933(E) HD 49933(H) 



[Na/H] 


+0.29 ± 0.04 


+0.06 ± 0.04 


-0. 


.14 + 0.04 


+0, 


.01 + 0, 


.03 


-0.02 + 0. 


,04 


-0. 


,32 + 0, 


.09 


-0. 


.35 + 0. 


.05 


[Si/H] 


+0.30 ± 0.03 


+0.05 ± 0.03 


-0. 


,07 ± 0.04 


+0, 


.05 + 0. 


.04 


+0.12 + 0. 


,05 


-0. 


,32 + 0, 


.04 


-0. 


.33 + 0. 


.04 


[S/H] 


+0.23 + 0.07 


-0.07 + 0.05 


-0. 


,12 + 0.03 


+0, 


.14 + 0, 


,20 


+0.02 + 0, 


,18 


-0, 


,25 + 0, 


,08 


-0, 


,31 + 0, 


,05 


[Ca/H] 


+0.23 + 0.03 


+0.02 + 0.03 


-0. 


,12 + 0.04 


-0, 


.04 + 0, 


,03 


-0.02 + 0, 


,06 


-0, 


,43 + 0, 


,04 


-0, 


,39 + 0, 


,05 


[Sc/H] 


+0.26 + 0.04 


+0.05 + 0.03 


-0. 


,13 + 0.08 


-0, 


.09 + 0, 


,14 


+0.14 + 0, 


,15 


-0, 


,35 + 0, 


,08 


-0, 


,34 + 0, 


,10 


[Ti/H] 


+0.26 ± 0.04 


+0.08 ± 0.05 


-0. 


,11 + 0.08 


-0, 


.12 + 0. 


,08 


-0.06 + 0, 


,04 


-0, 


,29 + 0, 


,04 


-0, 


,29 + 0, 


,04 


[V/H] 


+0.22 ± 0.04 


+0.07 ± 0.05 


-0. 


,21+0.09 


-0, 


.10 + 0. 


,07 


+0.03 + 0, 


,05 


-0, 


,43 + 0, 


,11 


-0, 


,50 + 0, 


,11 


[Cr/H] 


+0.27 + 0.03 


-0.02 + 0.04 


-0. 


,21 + 0.05 


-0, 


.09 + 0, 


,04 


-0.03 + 0, 


,03 


-0, 


,37 + 0, 


,05 


-0, 


,39 + 0, 


,08 


[Mn/H] 


+0.27 + 0.03 


+0.01 + 0.04 


-0. 


,25 + 0.05 


-0, 


,11 + 0, 


,08 


-0.08 + 0, 


,05 


-0, 


,29 + 0, 


,10 


-0, 


,43 + 0, 


,08 


[Fe/H] 


+0.24 + 0.02 


+0.02 + 0.02 


-0. 


,16 + 0.03 


-0, 


,02 + 0, 


,02 


-0.02 + 0, 


,02 


-0, 


,34 + 0, 


,03 


-0, 


,37 + 0, 


,03 


[Co/H] 


+0.19 ± 0.06 


-0.01 ± 0.05 


-0. 


,19 + 0.07 


-0, 


.02 + 0. 


,06 


+0.02 + 0, 


,07 


-0, 


,52 + 0, 


,05 


-0, 


,53 + 0, 


,07 


[Ni/H] 


+0.25 + 0.03 


+0.02 + 0.03 


-0. 


,23 + 0.04 


-0, 


.07 + 0. 


,03 


-0.04 + 0, 


,04 


-0, 


,53 + 0, 


,05 


-0, 


,51 + 0, 


,05 


[Cu/H] 


+0.24 + 0.06 


-0.02 + 0.05 


-0. 


,17 + 0.07 


-0, 


.08 + 0, 


,03 


-0.10 + 0, 


,08 


-0, 


,48 + 0, 


,07 


-0, 


,46 + 0, 


,08 


[Y/H] 


+0.31 ± 0.08 


+0.03 ± 0.03 


-0. 


,16 + 0.05 


+0, 


.00 + 0. 


,07 


+0.03 + 0, 


,04 


-0, 


,21 + 0, 


,08 


-0, 


,22 + 0, 


,07 


[Ba/H] 


+0.26 ± 0.06 


-0.05 + 0.07 


-0. 


,03 ± 0.08 


+0, 


.10 + 0. 


,11 


+0.23 + 0, 


,13 


-0, 


,17 + 0, 


,09 


-0, 


,20 ± 0, 


,08 


[Nd/H] 


+0.22 + 0.06 


-0.21 + 0.08 


-0. 


,08 + 0.06 


-0, 


.10 + 0, 


,12 


+0.00 + 0, 


,06 


-0, 


,25 + 0, 


,10 


-0, 


,21 + 0, 


,10 




HD 175726 


HD 177552 


HD 184663 


HD 46558 


HD 171834 


HD 49434 


HD 55057 


[Na/H] 


-0.07 + 0.06 


-0.10 + 0.10 


-0. 


,14 + 0.09 


-0, 


.06 + 0. 


,15 


-0.09 + 0, 


,16 


-0, 


,07 + 0, 


,18 


+0, 


,24 + 0, 


,19 


[Si/H] 


+0.01 + 0.04 


+0.15 + 0.10 


+0. 


,09 + 0.14 


+0, 


.15 + 0, 


,10 


+0.00 + 0, 


,20 


-0, 


,12 + 0, 


,15 


+0, 


,38 + 0, 


,25 


[S/H] 


+0.10 + 0.08 


+0.00 + 0.16 


+0. 


,21 + 0.07 


-0, 


.29 + 0. 


,18 


+0.01 + 0, 


,15 


+0, 


,11 + 0, 


,17 








[Ca/H] 


+0.06 + 0.04 


-0.19 + 0.09 


-0. 


,17 + 0.09 


-0, 


.05 + 0. 


,09 


-0.15 + 0, 


,08 


+0, 


,03 + 0, 


,10 


+0, 


,23 + 0, 


,17 


[Sc/H] 


+0.10 ± 0.07 






























[Ti/H] 


+0.07 ± 0.06 


+0.19 + 0.15 


-0. 


,04 + 0.09 


+0 


.02 ± 0, 


.14 


-0.06 ± 0, 


.10 


+0, 


.17 + 0, 


.22 


+0, 


.43 ± 0, 


.31 


[V/H] 


+0.12 + 0.08 






























[Cr/H] 


+0.03 + 0.04 


+0.00 + 0.04 


-0, 


,11 + 0.07 


+0 


.01+0, 


.09 


-0.03 + 0, 


,08 


+0, 


,07 + 0, 


,09 


+0, 


,27 + 0, 


.12 


[Mn/H] 


-0.01 + 0.04 


+0.10 + 0.12 


-0. 


,16 + 0.15 


-0, 


.07 + 0, 


,12 


+0.09 + 0, 


,18 














[Fe/H] 


+0.03 ± 0.03 


-0.03 ± 0.03 


-0. 


,13 + 0.06 


+0, 


.04 ± 0. 


.07 


+0.01 + 0, 


.07 


+0, 


.09 ± 0. 


.07 


+0. 


.26 ± 0. 


.08 


[Co/H] 


-0.02 ± 0.07 


-0.01 + 0.14 


-0. 


,10 + 0.13 


-0, 


.08 + 0. 


.18 


+0.02 + 0, 


.13 


+0, 


.13 + 0. 


.15 


+0. 


.20 + 0. 


.18 


[Ni/H] 


+0.01 + 0.04 


-0.08 ± 0.10 


-0. 


,19 + 0.14 


-0 


.02 ± 0. 


.17 


-0.02 ± 0, 


.11 


+0, 


.06 ± 0. 


.08 


+0. 


.38 ± 0. 


.11 


[Cu/H] 


-0.10 + 0.08 






























[Y/H] 


+0.16 ± 0.05 






























[Ba/H] 


+0.28 ± 0.08 


+0.16 ± 0.14 


-0. 


,26 + 0.13 
























[Nd/H] 


+0.20 ± 0.09 































